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Abstract 

We consider a Dirac field coupled minimally to the Mielke-Baekler 
model of gravity and investigate cosmological solutions in three di- 
mensions. We arrive at a family of solutions which exists even in the 
limit of vanishing cosmological constant. 
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1 Introduction 



Gravity in three dimensions has attracted a lot of attention since the inven- 
tion of the topologically massive gravity by Deser, Jackiw, Tempelton pQ-[3]. 
Another cornerstone is the discovery of the BTZ black hole solution to the 
Einstein's theory jl] . Since the inclusion of torsion makes the gravity models 
richer, Mielke and Baekler (MB) generalized the topologically massive gauge 
model of gravity by adding a new translational Chern-Simons term to the 
lagrangian of the standard topologically massive gravity [5] . There is a wide 
literature on the MB-model of gravity coupled to various material sources, 
see for example [5]- [ID] arid the references therein. We, however, realized that 
there is a tiny amount of work considering possible spinor couplings to the 
Einsteinan gravity and the standard topologically massive gravity in three 
dimensional spacetimes [H]-[T3]. For the MB-case the situation is worse 
[15]. Thus we intend to fill in this gap to gain new insights to cosmological 
problems in the context of relationship between spinor and gravity. 

We are interested in for investigating roles of a Dirac field in cosmology, 
because it is essential to a satisfactory description of relativistic fermions. 
Furthermore, a self-interacting or nonlinear Dirac field can yield negative 
pressure and thereby accelerate the early and the late-time expansion of the 
universe [16]. On the other hand, the authors show in the Ref.[T7] that 
a single Dirac field can give rise to inflation within the four dimensional 
Einstein-Cartan theory, and prove compatibility of the Dirac-field model with 
the observations by calculating the power spectrum of density fluctuations 
of the Dirac field. Because of the inflationary (exponential) nature of our 
solution fl33l) this work may be seen in favor of those results. 

The outline of the paper is as follows. Since we will be using the coordinate 
independent algebra of the exterior forms, in the Section [2] we fix our nota- 
tions and conventions for the Riemann-Cartan spacetimes and a Dirac spinor 
in three dimensions. In the Section [3] after introducing the topological gauge 
lagrangian of gravity and the Dirac lagrangian we obtain the field equations 
by varying independently the total lagrangian with respect to the orthonor- 
mal basis 1-forms, e a , the full connection 1-forms, u ab and the adjoint of the 
Dirac field, In the Section H] firstly we write down the homogenous and 
the isotropic metric in the plane polar coordinates and make torsion ansatz. 
Secondly we compute the Dirac equation and its adjoint. This allows us to 
calculate the Dirac energy-momentum 2-forms explicitly. Then we obtain all 
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field equations via the computer algebra system, the Reduce and the package 
Excalc [IB]) [IS]- Thus we determine a class of solution for the homogenous 
and isotropic geometry and a certain Dirac spinor. We summarize and dis- 
cuss the results in the Section |5j 

2 Notations and Conventions 

We specify the Riemann-Cartan space-time by a triplet (M, g, V) where M 
is a 3-dimensional differentiable manifold equipped with a metric tensor 

9 = Va b e a (8) e b (1) 

of signature (—,+,+). e a is an orthonormal co- frame dual to the frame 
vectors X a , that is e a (X ) = Lbe a = 5% where t& := tx b denotes the interior 
product. A metric compatible connection V can be specified in terms of 
connection 1-forms u a b satisfying Ub a = —u a b- Then the Cartan structure 
equations 

de a + u\ A e b = T a , (2) 
du a b + oo a c A uj c b = R a b (3) 

define the space-time torsion 2-forms T a and curvature 2-forms R a b, respec- 
tively Here d denotes the exterior derivative and A the wedge product. These 
tensor-valued 2-forms satisfy the Bianchi identities 

DT a = R a b A e b , DR a b = . (4) 

Because of the first Bianchi identity, the Ricci tensor does not need to be 
symmetric in the Riemann-Cartan spacetimes. We fix the orientation of 
space-time by choosing the volume 3-form *1 = e° A e 1 A e 3 where * is the 
Hodge star map. In three dimensional space-times with Lorentz signature 
for any p-form ** = —1. We will use the abbreviations e ab " := e a A e b A • • • 
and Lab-.- '■= L a Lb • • •- It is possible to decompose the connection 1-forms in a 
unique way as 

u\ = Gj\ + K\ (5) 
where u a b are the zero-torsion Levi-Civita connection 1-forms satisfying 

de a + u a b A e b = (6) 
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and K a b are the contortion 1-forms satisfying 

K\ A e b = T a . (7) 
The curvature 2-forms are also decomposed in a similar way: 

R a b = R\ + DK\ + K\ A K% (8) 

with 

DK\ = dK\ + Qj a c A K c b - Cj c b A K a c . 

We are using the formalism of Clifford algebra C£ 12 -valued exterior forms. 
C£\ t 2 algebra is generated by the relation among the orthonormal basis {70, 71, 72} 

7V + 7 b 7 a = 2r] ab . (9) 

One particular representation of the 7 a 's is given by the following Dirac 
matrices 

70 = ( -1 I ) ' 71 = ( I I ) ' 72 = ( I -1 ) ■ (10) 

In this Dirac spinor \I> can be represented by a 2-component column 

matrix. Thus we write explicitly the covariant exterior derivative of its 
Dirac conjugate and the curvature of the spinor bundle, respectively, 

D^> = cM + l -a ab ^uj ab , = - ^a ab u ab , D 2 ^> = l -R ab a ab ^ (11) 

where a ab := \[j a , 7&] = \^abc~f c are the generators of the Lorentz group 
50(1, 2). The Dirac adjoint is \l/ := \E^7o- We frequently make use of the 
identity 

IcOab + Cable = ^abc ■ (12) 

3 Field Equations 

The field equations of our model are obtained by varying the action 

I[e a 1 oo ab ^]= [ (L G + L D ) (13) 

J M 
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where the gravitational lagrangian density 3-form is given by Mielke and 
Baekler [S] 

L G = a -R ah A *e ab + A*l + h -T a A e a + °-{uj a b A du\ + 2 -uj\ A u\ A u; c a ) (14) 
and the (hermitian) Dirac lagrangian density 3-form 

L D = ^(f*7AM-MA *7*) +im***l. (15) 

Here the gravitational constants a, mass m and the Dirac field \l/ have the 
dimension of length -1 , the gravitational constant b has the dimension of 
length -2 , c is dimensionless constant, and the cosmological constant A has the 
dimension of length -3 . The case b = corresponds the topologically massive 
gravity with cosmological constant. The hermiticity of the lagrangian f JT5|) 
leads to a charge current which admits the usual probabilistic interpretation. 
|T a A e a is known as the torsional Chern-Simons term which corresponds the 
usual Chern-Simons 3-form for the curvature, (l/2)(u a b A du b a + (2/3)u a b A 
u b c A u c a ), for the curvature. 

We obtain the field equations via independent variations with respect to 
e a ,u ab , Thus e a -variation yields the FIRST equation 

-^e abc R bc -\*e a -bT a = r a , (16) 
w a6 -variation yields the SECOND equation 

- a -e abc T< + b -e ab + ciU = S ab , (17) 
and ^/-variation yields the Dirac equation 

*7A(D-^V)* + m**l = 0, (18) 

where T, ab = — j^^e ab is the Dirac spin angular momentum 2-form and r a 
is the Dirac energy-momentum 2-form 

r a = l -*e ba A [* 7 6 (Dtt) - (Dtt) 7 6 tt] + imM*e fl . (19) 
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For future convenience by using the Dirac equation (TTS)) and its conjugate 
(D — \V)^ A *7 — m\I/*l = we rewrite the Dirac energy- momentum and 
the spin angular momentum 2-forms as 

r a = ~\*jb(d a *)-(d a W)>y b V]*e b + Scj beta e be (20) 
S afe = -Se ab (21) 

where d a := i a d, u bc , a := L a u bc and S := |\I>\I/. 



4 A Class of Cosmological Solution 

We consider the homogeneous and isotropic metric in the plane polar coor- 
dinates (t, r, 0) 

g = -dt 2 + ^^-Ar 2 + r 2 Q 2 {tW (22) 
1 — kr z 

with the expansion factor Q(t) and curvature index k. We choose the or- 
thonormal basis 1-forms 

e° = dt , e 1 = -B^=dr , e 2 = rQ(t)d0, (23) 

leading to the Levi-Civita connection 1-forms 

-o Q i -o Q 2 -l \f\--kr 2 2 fr,A\ 

u\ = —e , w u 2 = — e , u 2 = e (24) 

Q Q rQ 

where dot denotes the derivative with respect to t. Two independent func- 
tions are enough to describe the most general torsion preserving homogeneity 
and isotropy of the spacetime [20] . Thus we choose the torsion 



T° = u(t)e m + v(t)e 12 , T 1 = v(t)e 02 , T 2 = -v(t)e 01 - u(t)e 12 . (25) 
We then calculate the contortion 1-forms via (jJJ) 

A' 1 = ne°-V, ^° 2 = V, ^ 2 = - M e 2 + V. (26) 
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Now we write explicitly the full connection 1-forms with the help of the 
equation (j5J) 



o Q 1 , v 2 
u i = -ue - —e + -e 



v i Q 2 
c 02 = - 



W12 = -e - | 



(27) 



Thus the curvature 2-forms are computed as 
(4Q + t? 2 Q)e 01 - 2i>Qe 02 + 2wQe 12 



4Q 

2rvQe 01 - Au^/Y^W 1 + r(4w 2 - t> 2 )<3 - 4rQ e 02 - 4rwQe 



(2* 



4rQ 



2ruvQ 2 e m - AruQ 2 e 02 + 4rQ 2 + rw 2 Q 2 + Akr - 4uQV~L - kP 



,12 



4rQ 2 



For the consistency of the equations, that is, for keeping all Q, u, v functions 
to be dependent only on the cosmic time t we have to assume \P depending 
on both t and r such that \& = r) = £{t)/y/r. Thus the Dirac equation 
( fl8|) and its conjugate turn out to be 



£ = -Q/Q + (3t;/4 - m) 7o j £ 
t = t\-Q/Q-(3v/4-m) l0 



(29) 
(30) 



These results enable us to determine the Dirac energy-momentum 2-forms 
(EOD 

t = 2S\-ue m + (2m - v)e 12 ) , 
Tl = -l[2Qe 01 + vQe 02 }, 
S 



r 2 = ^[rvQe 01 _ 2rQe 02 - 2{ruQ + Vl - kr 2 )e 12 ] . (31) 
Consequently we obtain the following class of solution to the field equations 

(USD and (ED 

ab — 2cA 



u — , u 



a 2 - 26c 



5 = 



(32) 
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Q(t) 



1 

2h 



(33) 



where ci is a constant and /i 2 = [4(a 2 — 2bc)(b 2 — aX) — (ab — 2cA) 2 ]/4(a 2 — 
26c) 2 . If k < and ci = (In \k\)/(2h), the space-time characterized by 
(|33|) which recalls the inflationary solution of the General Relativity in four 
dimensions has a singularity, i.e. Q(0) = 0. Here we also notice that it 
must be a 2 — 2bc ^ for the solution, but even if A = we have a solution. 
The condition a 2 — 2bc ^ was shown to be arising from the demands of 
a canonical constrained hamiltonian analysis in [2T]. For more readings one 
can consult the equation (38) and the succeeding remarks of that work with 
the replacement a — > a, — > 2c and — > b. However, a prescription to 
approach the singular point a 2 = 2bc and the canonical structure analysis of 
the theory in the space of parameters is discussed in 



(34) 



The next task is to solve the Dirac equation ([29 
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-Q/Q 3v/A - m 
-3v/4 + m -Q/Q 



6 
6 



We can decouple these equations by virtue of a unitary transformation 

£ = Up with U = -4= ( 1 1 
Thus we arrive at the set of decoupled equations 



(35) 



Pi 



-Q/Q + i{3v/4- m 
The solutions can be written as 
Pi(t) 



Pi > P2 



-Q/Q-i(3v/4-rn) p 2 . (36) 



p+ e +m) 



Q(t) 



P2(t) 



Q(t) 



-i<t>(t) 



(37) 



where p± are the integration constants and <fi = (3^/4 — m)t is the phase 
factor. Then we calculate the Dirac condensate that is zero according to 



ip+i 



4rQ 2 



0. 



(38) 



Consequently it must be |p+| 2 = |p~| 2 - Furthermore since S = 0, the Dirac 
spin 3-form is zero, = 0, via (|2~T1) as well. 
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5 Concluding Remarks 



In this work we considered the minimal coupling of a Dirac particle to the 
MB-model of gravity in three dimensional Riemann-Cartan spacetime. After 
computing the variational field equations we considered the homogenous and 
isotropic metric, and made an ansatz for torsion. Then we calculated the 
Dirac equation which gave us the opportunity of obtaining the Dirac energy- 
momentum 2-forms. Consequently we were at a position of writing all field 
equations explicitly. Thus we obtained a family of solution f )32|) . fl33|) and 
fl37|) which is valid for a 2 — 26c 7^ and even if A = 0. 

It is worthwhile to pay more attention to the result, S = 0, given by the equa- 
tion (132 p . It may imply that even though the Dirac field has no effect on the 
geometry, the later affects the time evolution of the spinor field. Thus our 
solution may evoke the Einsteinian cosmological solution satisfying perfect 
cosmological principle in which A 7^ is a crucial condition. However, our 
model and the geometry are totally different from Einstein's theory and the 
Riemann spacetime, respectively. In fact, since the presence of a Dirac spinor 
modifies the geometry of space-time by affecting the connection that deter- 
mines the notion of parallel transport, it is no longer possible to conclude 
that the geometry is totally unaware of the presence of a Dirac spinor. 

Meanwhile, there are works in which the Einstein-Cartan-Dirac (ECD) the- 
ory was investigated in four dimensions, for example [22] , [21] , [22] • Even 
they were done in four dimensions, they are related to our work, especially, 
with regards to S = 0. For clarifying the relation we firstly notice that 
5 = means *(S ab A *S afe ) = equivalent to S abc S afec = in the language of 
tensors with the notation S a b := S a fc c *e c . In four dimensions because of the 
complexity of the ECD problem the authors of those papers proceed with 
the following strategy. They firstly determine orthonormal coframe e a by 
solving the Einstein equation (with cosmological constant) in vacua. Then 
they give the constraint by hand E a & c £ abc = for Finally they compute 
torsion and spinor through the field equations. For massless Dirac spinors 
this constraint corresponds to the helicity state equation |(1 ± ^75)^ = ^ 
[25] . The neutrinos are considered as massless fermions and correspondingly 
have definite helicities according to the standard model. Therefore, they may 
be treated in this context. Thus in the Ref.[26] even though the authors did 
not follow directly the strategy outlined above, their results and calculations 
are parallel to it. However, in three dimensions because of relative simplicity 
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one does not need to assume any constraint like that from the beginning. We 
encounter with that as a result of the equations. But this is not a generic 
result, that is, for some other models and their solutions in (l+2)-dimensions 
it may be S afec S afcc « S 2 ^ 0, see e.g. [H]. 



Appendix: Irreducible Decompositions 

In this section we give briefly the irreducible pieces of torsion, contortion and 
curvature in three dimensions. Firstly torsion which has nine components can 
be decomposed 

rpa _ (I) rpa _|_ (2)rpa _|_ (3) rpa j-gg-j 
#9 #5 #3 #1 

where < 2) T a = -\{i b T b ) A e a , ^T a = \i a (e b A T h ) and ^T a = T a - ^T a - 
(3)j>a_ j n ^-^jg sec tion the notation with the number under a brace is for the 
number of components of that part. They have the properties, ^T a A e a = 
(2) T a A e a = and L a {1) T a = t a (3) T a = 0. Our choice p)} corresponds to 

/ e 01 \ / e 12 \ 

(i) T a = Q) (2) T a = / o , ^T a = v e 02 (40) 

which means 4 = 3 © 1. After the solution ( 1321) we are left only with ( 3 )T a . 
Secondly one can decompose the contortion having nine components 

^ = ^^ + ^U + ( ^JU ( 41 ) 
#9 #5 #3 #1 

where ^K ab = \[e a A (6 c K cfe ) - e fe A (^ c K ca )], (3) if ab = -\i ab {K cd A e cd ) and 
WiCrf, = K ab -^K ab -^K ab . They have the properties t a ^K ab = t a ^K ab = 
and ^K ab A e afc = ^K ab A e ab = 0. For our case fl26|) we possess again 
(1) K ab = 0, but nonzero ^K ab and ^K ab , i.e. 4 = 3 © 1. Besides after the 
solution ([3"2l only ^K ab survives. 

Finally one can split the curvature with nine components 

^= { ^B^ + ^R^ + ^R^ (42) 

#9 #5 #3 #1 



9 



where {2) R ab = \(e a A t b - e b A L a )(e c A R c ), {3) R ab = \Re ab and {1) R ab = 
Rah — ^Rab — ^Rab with R a = i b ' R ba and R = L a R a . They have the properties 
i ab WR ab = i ab {2) R ab = 0, WR ab A e b = ®R ab A e b = and e b A t a ( 1 ^ a6 = 0. 
For the solution fl32l) although all three pieces of curvature and of torsion are 
nonzero, since DT a = the Ricci tensor is symmetric. 
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